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FACTORS INFLIJKNCING THE GUIDANCE ACCURACY 


OF A PLANETARY PROBE 
By Thomas 3. Murtagh 


SUMMARY 


The influence of the type (i.e., relative range, range rate, or 
both) and accuracy of spacecraft onboard radar data, nominal entry param- 
eters for the probe trajectory, and planet radius uncertainty for un- 
manned planetary probe guidance is presented. The probe is assumed to 
be deployed from a manned spacecraft at the destination planet's rphere 
of influence, with the inclination of the probe trajectory equal -o that 
of the spacecraft approach hyperbola. The results of the analysis indicate 
that a probe entry guidance corridor of the order of 20 n. mi. can be 
obtained for a total midcourse AV less than 100 fps , with spacecraft 
onboard radar range data contributing more information to entry corridor 
determination than range-rate data. These results were more sensitive 
to nominal entry trajectory parameter variations than to variations in 
the assumed error in the radius of the target planet. 


INTRODUCTION 


The analysis presented in this note extends and investigates some 
of the assumptions of reference 1 in more detail. The most significant 
assumption was that the radar onboard the manned spacecraft could measure 
precisely the type of data - relative range and/or range-rate - to the un- 
manned probe. That is, the onboard estimates of the state vectors for the 
probe and spacecraft have equivalent accuracy. The nominal probe entry 
trajectory was also specified so that the probe would arrive at the 
nominal entry altitude approximately one hour before the spacecraft 
arrived at periapsis of the approach hyperbola. These nominal entry 
parameters were fixed in order to simplify the analysis. However, the 
entry parameters are a function of the type of probe mission considered. 
For example, shallow entry flight-path angles would be desired for probe 
missions which require a "skip-out" of the atmosphere into an orbit about 
the planet; impact-type probes require much steeper flight-path angles. 

For the analysis presented, a reference probe mission was chosen for 
determining the effect of the type and accuracy of spacecraft onboard 
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radar data on the guidance corridor. Variations in planet radius error 
and nominal entry trajectory parameters are presented in order to illustrate 
their effect on the entry corridor. 

The results are presented for a conjunction class Mars mission and 
can be considered valid for other missions with similar characteristics. 
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SYMBOLS 


sensitivity vector which relates star /planet-horizon 
included angle deviations to spacecraft state vector 
deviations 

sensitivity matrix which relates spacecraft to piobe 
relative range and range-rate deviations to spacecraft 
to probe state vector deviations 

6x3 compat ability matrix defined by 
probe uncertainty covariance matrix 
spacecraft uncertainty covariance matrix 
6x6 giii dance matrix 

3 3 sub-matrices of the guidance matrix 

3 X 12 sensitivity matrix defined by equation (9) 

identity matrix of appropriate dimensions 

weighting matrix defined in equation ( 13 ) 

3x3 matrix defined by equation (16) 

3x3 matrix defined in equation (l 3 ) 

covariance matrix of velocity correction execution 
error 

12 X 12 uncertainty covariance matrix for coupled 
spacecraft /probe system 
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covariance matrix of measurement errors defined by 
equation (lU) 

planet radius 

probe position vector with respect to planet 
spacecraft position vector with respect to planet 
magnitude of spacecraft position vector 
unit vector to a star 
current time 

velocity vector of probe with respect to planet 

velocity vector of spacecraft with respect to pla.iet 

spacecraft /probe relative velocity vector (* V - V ) 

P 8 

probe dispersion covariance matrix 

spacecraft dispersion covt -iarice matrix 

star /planet-horizon included angle 
6x6 probe state transition metrix 
small deviation from reference value of ( ) 
degradation to ^^(t) as result of spacecraft /probe 
reparation maneuver 

12-dimensional augjnented state vector 

one-half planet disc subtended angle (sin 6 = r_/r ) 

B s 

12 X 12 augmented state transition matrix defined by 
equation (3) 


u 



p(t) 

^(t) 



4>(t ,t) 


relative position vector of probe with respect to 

spacev'raft (p ■ r - r ; 

p s 

magnitude of relative position vector 
spacecraft /probe relative range-rate ^ 

variance of star /planet-horizon observable 

onboard radar range variance 

onboard radar range-rate variance 

projected time 

6x6 spacecraft state transition matrix 


Si’^-ercripts : 

( y ( ) after navigation measurement or guidance maneuver 

( )~ ( ) before navigation measurement or guidance maneuver 

( )" Inverse of ( ) 

( ) Transpose of ( ) 

ANALYSIS 


Assiimptlons 


/' * 


Following are the postulates for this analysis: 

1, The unmanned probe is deployed from the manned spacecraft at 
the destination planet sphere of influence (SOI) so that the inclination 
of the probe trajectory equals that of the spacecraft approach hyperbola. 

2. Spacecraft and probe position and velocity uncertainties are re- 
duced by simultaneously processing onboard radar relative range and/or 
range-rate data and star /planet -horizon included angle data from the 
spacecraft using a Kalman filter. 
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3. Conic reference trajectories were assumed, and the state transi- 
tion matrix used to propagate the errors was derived analytically for 
twc-tody conic trajectories. 

4. Variable time of arrival (VTA) guidimce logic was used to compute 
the root-mean-square (RMS) velocity corrections for the probe and the 
corresponding vacuum periapsis radius dispersions. The probe entry 
guidance corridor is computed by multiplying the radius dispersion by a 
factor of six (± 3a about the nominal). 

5. A summary of the navigation and guidance system nominal errors 

is presented in table I. The spacecraft onboard radar errors are consist- 
ent with the values used in reference 2. 


Reference Mission 

The reference mission chosen for the eincdysis presented in this 
note was a 1977 Mars stopover (ref. 3) which has an outbound trip time 
of 360 days, a Mars orbit time of 300 days, and a return to Earth time 
of 320 days. The reference trajectory characteristics for this mission 
are summarized in table II. 

The Earth-injection covai-iance matrix for this trajectory was 
diagonal with RMS position and velocity errors of U n. mi. and 16 fps. 
Earth-based radar range and range-rate data were processed during the 
departure phase of the mission, i.e., within the Earth SOI, and was 
followed by onboard sextant planet-star tracking in the heliocentric 
phase, i.e., between the Earth and Mars SOI’s. For these onboard 
raeastirements the sighting bod^' was selected according to the optimality 
criteria outlined in reference U. The error model for the Earth-based 
radar system is discussed in reference 5; the onboard sextant error 
model can be found in reference 1. 

The Earth-bused radar and onboard optical navigation data were 
processed every 30 minutes and each half-day during the departure and 
heliocentric phaises , respectively. Three midcourse fixed time of arrival 
(FTA) velocity corrections requiring a total AV of 7^.12 fps were ex- 
ecuted to update the spacecraft dispersion matrix prior to the unmanned 
probe deployment at the Mars SOI (approximately 312 000 n. mi. from the 
planet ) . 


Spacecraft /Probe Navigation and Guidance System 

The spacecraft /probe tracking geometry is illustrated in figure 1. 
For the study presented in this note it was assumed that the spacecraft 
onboard radar could measure the relative range and/or range-rate to the 
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probe and simultaneously use an onboard optical sensor, l.e., sextant, 
to measure the Included angle between the Mars horizon and a star. This 
procedure seems feasible since the onboard radar can track the probe 
continuously and, when the spacecraft horizon-star measurement Is fed 
Into the onboard computer, a command could automatically be set up In 
the navigation program which would call for simultaneous data processing 
of the radar range and/or range-rate information. 

The navigation data can be processed In the onboard computer vulng 
a Kalman filter. The structure of the filter equations for the coupled 
spacecraft /probe system is identical to the conventional Kalman equations 
but has increased state vector dimensions (ref. 2 and 6). For the prob- 
lem considered here, the state vector is 12-dlmensional and includes the 
spacecraft position and velocity as well as the position and velocity of 
the unmanned probe. The equation which relates deviations in this state 
vector at time, t, to deviations at time, t, is 
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then equation (l) becomes 


6^(1) * 0(t ,t ) 6?(t ) (U) 


The initial 12 x 12 covariance matrix for the coupled system, i.e., at 
spacecraft /probe separation, is defined by 


P(t ) 
o 




E„(t ) 
P o 


(5) 


where is spacecraft uncertainty covariance matrix and 


E^(t^) *= E^(t^) + AE(t ) (6) 

p o so o 

The term AE(t ) is the degradation to the probe uncertainty matrix as a 
result of the^assumed imperfect separation maneuver. The equation for 
propagating this matrix between measurements is given by the expression 

P(t) = 0(T,t)P(t)0'^(T,t) ( 7 ; 


The equation which relates deviations in the observables, i.e., star 
horizon included angle, range, and range-rate, to state vector deviations 
is 


6B(t) 

6p(t) 

6p(t) 


H(x)6t(T) 


^8) 


where the 3 ^ 12 matrix H(t) is written in partitioned x'orm as 


A(t ) 0 

-B(t) B(t) 


H(t) 


(9) 


8 


The vector A(t) Is defined by 


A(t) - 




(10) 


and the 2x6 matrix B(t) is given by 


B(t) = 


ap 

ap 




( 11 ) 


The partial derivatives required in equations (lO) and (ll) can be com- 
puted from the following relationships (fig. l). 
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(12c) 
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The equations required to update the uncertainty matrix, P(t), at the 
time of a measurement can now be written 


>*(t) = 


I - K(t)H(t) P"(t) 


] 


K(t) = P"(t)H^(t)m“1(t) 


M(t) = H(t)p“(t)H^t)+R(t) 




( 13 ) 


where the 3 x 3 covariance matrix of measurement errors, R(t), is defined 
by 


R(t) = 




ilh) 


and the (-) and (+) superscripts refer to a quantity before and after the 
measurement (or correction) respectively. 

If it is assumed that the navigation and guidance systems are un- 
coupled, then the spacecraft and probe dispersion matrices are propagated 
separately using the equations 


Xg(i) = $(T,t)X^(t)4>'^(T,t) 




X (t) = r(i,t)x (t)rMi,t) 

It h' 


( 15 ) 
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When reeisonable confidence is obtained in the trajectory estimates of 
either the probe or spacecraft, guidance maneuvers are commanded for the 
appropriate vehicle to restore the dispersed trajectory to specified 
nominal conditions. For example, if a guidance correction is commanded 
for the probe at time, t, then the RMS estimate of the required AV is 
computed from the square root of the trace of the equation 

L(t) . G2(t)] |^Xp(T)-P^(T)j |g^(t) 


where G^(t) and are subraatrices of the guidance matrix G(t) discussed 

in references 1 and 7. The matrix P^(t) is a submatrix of the augmented 
uncertainty matrix P(x) defined by 


P(t) 


P^(t) 

PgCt) 

P3(t) 

P^(t) 


(17) 


after one or more navigation measurements are processed. 

The probe uncertainty and dispersion matrices are modified according 
to the equations (ref. l) 


= Pij'(t) + CN(t)c''' 

*(t) = jl+G(T)j j^Xp*(T)-P^‘(T) j |i+G(t) j'^+P^*(T) ) (18) 




Equations similar to (l6) and (l8) are used to calculate the RMS AV and 
matrix updates for the spacecraft if ^^(t) and P|^(t) are replaced by 

X (t) and P_(t) respectively. The covariance matrix of velocity correction 
s ^ 

execution error, N(t), is derived and discussed in reference 8. 


RESULTS 


Nominal Probe Trajectory 

Assuming that the probe is deployed from the spacecraft at the Mars 
SOI so that the inclination of the probe trajectory equals that of the 


• /' 
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spacecraft approach hyperbola, three entry parameters remain to be 
specified - altitude, speed, and flight-path angle. The altitude 

assumed for the probe nominal trajectory was 315 000 ft. Figure 2 presents 
plots of probe separation AV versus entry speed for entry flight-path 
angles of 0° and -U5° and for this entry altitude. The minimum separation 
velocity occurs when the angle between the separation AV vector and the 
spacecraft velocity vector is 90°. From figure 2(a) it can be seen that 
the effect of variations in entry flight-path angle becomes more pronounced 
as the minimum of the curve is approached and decreases quite rapidly on 
either side of the minimum. The time required for the probe to reach 
vacuum periapsis (not shown on the figure) is inversely proportional to 
the entry speed. For the nominal probe trajectory an entry speed of 
l8 350 fps was chosen, with an entry flight-path angle of -5°* These 
values of entry parameters cause the probe to reach its vacuum periapsis 
approximately 20 minutes ahead of the spacecraft arrival at the periapsis 
of the approach hyperbola. The separation AV required was U5 fps. 

The spacecraft /probe relative range is plotted against time from 
separation for the nominal probe trajectory in figure 3. This relative 
range as a function of time depends upon the separation AV which, in 
turn, is a function of the specified entry parameters. The nominal probe 
trajectory chosen is very dependent on the relative range since the space- 
craft onboard radar must have a maximum range beyond which no tracking 
is possible. The maximum range to the probe for the above nominal tra- 
jectory was about 2800 n. mi. so that effective probe tracking could be 
assumed throughout the probe delivery phase of the mission, if the onboard 
radar range capability was at least that value. 


Influence of the Type of Onboard Radar Data 


The effect of the type of spacecraft onboard radar data and accuracy 
on the unmanned probe navigation is illustrated in figure U. In figure 
U(a) it is assumed that the onboard radar processes both range and range- 
rate data every 30 minutes for three sets of range and range-rate errors. 
The profile of the curves in this figure is not entirely what one would 
expect. The apparent anomaly in the data occurs between 10 and 50 hours 
from separation. Within this time span the larger radar errors produce 
lower projected vacuum periapsis radivn uncertainties than the corre- 
sponding smaller radar errors . The explanation for this phenomena is 
related to the correlation which exists in the uncertainty matrix for the 
coupled spacecraft /probe system. It should be pointed out that the data 
weights, i.e., K(t), computed by the filter are a function of both the 
spacecraft and the probe imcertainties propagated from a previous meas- 
urement. In the region 10 - Uo hours from separation, the a = 200 ft and 


o* = 2 fps radar errors produce lower projected probe uncertainties than 
= 100 ft and = 1 fps radar errors. The spacecraft errors, however, 
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projected to the periapsis of the approach hyperbola (data not shown) 
during this same time interval are slightly larger for the higher radar 
errors, i.e. , = 200 ft and * 2 fps , compared to values for the 

0 * 100 ft and o* * 1 fps radar errors. This implies that the filter, 

P P 

in this time interval, gave more weight to the probe data compared to the 

simultaneously processed spacecraft data. A similar effect occurs between 

Uo and 50 hours from separation for the o = 300 ft and o» « 3 fps radar 

P P 

errors. The explanation is equivalent to that presented above. 

Figures U(b), U(c), and U(d) present probe navigation data comparing 
the three possible radar data type combinations , again assuming that 
data was processed every 30 minutes . In these figures the range and 
range-rate curve was generated by simultaneously processing range and 
range-rate information in the filter. The range curve was calculated 
by processing only range data in the filter in a separate simulation 
run. A third computer run was required to generate the range-rate curve 
by processing only range-rate data in the filter described by equations 
(13). 


In figure U(b), between 28 and 57 hours from probe deployment, the 
range-only tracking provides a lower vacuum periapsis radius uncertainty 
than the combination of range and range-rate tracking. This effect is 
due to the low information content of the range-rate data which tends to 
degrade the range and range-rate combination. This effect is present to 
a lesser degree in figures U(c) and U(d). A possible explanation is that 
the relative weight of the range-rate data compared to the range data 
has diminished as the errors were increased. 

The probe guidance plots associated with the navigation data in 
figure U are presented in figure 5* Figure 5(a) illustrates the effect 
of the radar errors on the probe entry guidance corridor when range and 
range-rate data are processed every 30 minutes. Assume for the moment that 
a vacuum periapsis radius dispersion of U n. mi. is desired (corridor == 

2k n. mi.). The nominal radar errors achieve this specified dispersion 
for an RMS midcourse AV equal to 80 fps. Increasing the radar errors 
by a factor of two approximately doubles the AV required to achieve this 
same dispersion. Increasing the radar errors by a factor of three, however, 
does not triple the required AV. 

Figures 5(b), 5(c), and 5(d) illustrate the effect of the type of 
radar tracking on the guidance corridor for the three sets of radar 
errors. In these plots it is obvious that range-rate data alone will 
not allow the probe to achieve the specified 2U-n. mi. guidance corridor. 

An examination of figure 3 verifies this result. In that figure it is 
evident that the rate-of-change of the relative range is constant be- 
tween the time of probe deployment and 55 hours from separation, implying 


13 


an insensitivity of the range-rate measurement to probe state vector 
variations. And, in all the csises shown, the range-only tracking 
produces this corridor for less AV than the range and range-rate tracking 
combination. 


Influence of Mars Radius Uncertainty 

The effect of the Mars radius error on the probe entry guidance 
corridor is illustrated in figure 6. Nominal radar errors are assumed, 
with range and range-rate measurementr processed simultaneously every 
30 minutes, using the nominal probe reference trajectory previously dis- 
cussed. If a guidance corridor of 2U n. mi. is desired (vacuum periapsis 
radius dispersion = U n. mi.) then the midcourse AV required are 30, 60, 
and 100 fps for Mars radius errors of 2, 10, and 20 n. mi. 


Influence of Varying Nominal Probe 
Trajectory Parameters 

The effect of nominal trajectory parameter variation on probe 
guidance accuracy using nominal Mars radius and spacecraft onboard radar 
errors is presented in figure 7. Navigation data was processed every 
30 minutes . 

In figure T(a) the probe midcourse AV plotted as a function of 
the RMS vacuum periapsis radius dispersion for three values of entry 
altitude. The other entry parameters were fixed at their nominal values. 
As the entry altitude increases, the separation AV increases and the 
time from separation to vacuum periapsis decreases. The midcourse AV, 
for a given vacuum periapsis radius error, also increases as the entry 
altitude increases. It should be pointed out that the maximum relative 
range between the spacecraft and probe also increases with increases in 
the entry altitude. Consequently these guidance results are valid only 
if the onboard radar range capability is equal to or greater than 
7500 n. mi. Otherwise effective probe tracking cannot be assumed through 
out its delivery to the entry interface. 

Variations in the entry flight-path angle are considered in figure 
7(b). As the flight-path angle decreases, the separation AV increases 
slightly and the time to the probe vacuum periapsis remains almost con- 
stant. For a specified corridor the midcourse AV decreases as the 
flight-path angle is varied from -5*^ to - U0°. 

Finally, variations in entry speed are illustrated in figure 7(c). 
Decreasing the entry speed produces corresponding reductions in the mid- 
course AV necessary to achieve a specified corridor. 


lU 

Spacecraft Navigation and Guidance 

All of the data discussed in the preceding sections involved the 
probe navigation and guidance. The spacecraft from which the probe is 
deployed must also be maneuvered to some safe target dispersion. How- 
ever, in order to keep the dimensions of this note within reasonable limits 
only a brief sketch of spacecraft data will be presented (fig. 8). 

Figure 8(a) presents the spacecraft navigation data from space- 
craft/probe separation to the approximate time of spacecraft arrival at 
the periapsis of the approach hyperbola. Similar data was presented in 
reference 1. The solid curve represents the projected periapsis radius 
uncertainty assuming that no probe tracking occurs , with Mars-horizon/star 
included angle measurements processed eveiy 3C minutes. The dashed curve 
presents analogous data assuming that the spacecraft now tracks the probe. 
It should be noted that better navigation results are obtained when the 
spacecraft tracks the probe than when it does not. This is a result of 
coupling the spacecraft /probe system through the measurements of relative 
range and range-rate. ITie effect noted here is quite analogous to re- 
ducing landmark location error in orbital navigation problems. There a 
spacecraft in orbit about a planet makes optical sightings to poorly 
defined landmarks and the output of the data processing system reduces 
the landmark position uncertainty els well as the uncertainty associated 
with the state of the orbiting spacecraft. 

The spacecraft gtiidance results are presented in figure 8(b). 

Again an improvement is noted when the spacecraft tracks the probe as 
compared to the case when it does not. For example, if a la periapsis 
radius dispersion of 3 n. mi. is specified, a midcourse AV of 13 fi)s is 
required if the spacecraft is tracking the probe compared to 25 fps if 
it is not tracking. 

The lower values of spacecraft AV for a specified periapsis radius 
dispersion compared to the probe results presented in the previous sec- 
tions is due to the fact that the initial spacecraft errors are smaller 
than the probe errors. [See equation (6).] If there were no probe de- 
ployment execution errors, i.e., AE(t^) = 0, then the spacecraft and 
probe guidance results would be approximately the same. 


CONCLUDING REMARKS 


The influence of the type and accuracy of spacecraft onboard radar 
data, nominal entry trajectory parameters, and planet radius error for 
unmanned planetary probe guidance has been presented. The results of 
the study indicate that a probe guidance corridor of about 20 n. mi. 


can be obtained for a total mldcourse AV around 100 fps , and that the 
spacecraft onboard radar relative range data contribute more information 
to entry corridor determination that the relative range-rate data. For 
the cases run, the corridor results were more sensitive to nominal entry 
trajectory parameter variations than to variations in the assumed un- 
certainty in the planet r^aius . 
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TABLE I.- NOMINAL 1 o PMS GiN ERRORS 
(a) Navigation system 


Onboard radar 

Range, ft 

Range rate, fps 

Onboard sextant , arc sec . 

Mars radius, n. mi. . • • 


100 

1 

10 

10 


(b) Guidance system 

Proportional percent 

Pointing, deg 

Cutoff, fps 


1 

1 


0.5 
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TABLE II.- CHARACTERISTICS OF THE 1977 MARS 
STOPOVER REFERENCE MISSION 

Launch from eeu*th, Julian date 2 UU3 UUO 

Earth injection velocity magnitude, fps 12 652 

Outbound trip time , days 360 

Mars stopover time, days 300 

Return trip time, days 320 

Periapsis altitude at Mars, n. mi.^ 200 

Entry velocity at eeu-th, fps 38 U63 

a 

The Mars approach hyperbola periapsis is located in the 
Martian Southern hemisphere. 
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Figure 1.- Spacecraft/probe tracking geometry. 
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Figure 3.- Spacecraft/probe relative range as a function of time from separation. 
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Figure A,- Unmanned probe navigation. 
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(b) Effect of onboard radar data type, CTp= 100 ft. 




(c) Effect of onboard radar data type, o« = 200 ft, ojp - 2 .0 fps. 
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(d) Effect of onboard radar data type, Op= 300 ft, - 3 .0 fps. 
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(c) Effect of onboard radar data type, <Tp = 200 ft, 2 .0 fps. 
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Figure 6,- Effect of Mars r»Hius error on probe guidance. 



Figure 7.- Influence of nominal trajectory parameter variation on probe guidance. 



Entry altitude 315 000 ft 

Entry flight-path angle, deg Entry speed 18 350 fps 


52 



sdj 'AV asjnoopjui aqojd glAld 


Entry altitude 315 000 ft 

Entry flight-path angle "5 deg 


* 


• • ✓ 


55 



soj 'AV ssJncDpiui aqojd gwa 


No probe tracking 


34 



Figure 8.- Spacecraft errors. 
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